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Abstract. MACRO can detect three topologies of neutrino 
induced events, corresponding to different parent neutrino ener- 
gies. The most numerous sample is made of upward throughgo- 
ing muons induced by atmospheric neutrinos of average energy 
100 GeV. The ratio of the observed over the expected events 
is 0.74 ± 0.036stat ± 0.046^^5 ± O.lStheor- The observed zenith 
distribution does not fit the expected one in the no-oscillation 
hypothesis, giving a maximum probability of 0.1%. Con- 
sidering the — Vt oscillation hypothesis the best probability 
(17%) is obtained for maximum mixing and Am^ of a few times 
10"^ eV^. The other detected samples are due to internally pro- 
duced events and upward-going stopping muons, corresponding 
to an average neutrino energy of around 4 GeV. The results con- 
cerning these lower energy sample show a deficit and a shape of 
the angular distributions in agreement with that predicted by 
the oscillation model suggested by the higher energy sample. 
No evidence of a neutrino signal due to dark matter particles in 
the direction of the core of the Earth and of the Sun has been 
found among the background due to atmospheric neutrinos and 
limits have been set. The neutralino hypothesis is investigated 
and limits on its mass are given. 



1. Introduction 

In the past years, the water Cherenkov experiments 1MB and Kamiokande 
have observed around the 60% of the expected double ratio of contained 

^ E-mail: montaruli@ba.infn.it 

1 



muon neutrino to electron neutrino interactions More recently, the 
Soudan 2 experiment has confirmed the anomaly using an iron-based de- 
tector On the other hand, the earlier results from the iron calorimeters 
NUSEX and Frejus are consistent with the expected number of contained 
events though with smaller statistics [||. 

Here we present the measurement of the high energy muon neutrino 
flux and its interpretation in terms of neutrino oscillations and the first 
results concerning the low energy events in MACRO. 



2. Neutrino detection with MACRO 

The MACRO detector important features for neutrino detection are its 
large area, fine tracking granularity, symmetric electronics with respect to 
upgoing versus downgoing particles, good timing and angular resolutions. 
Moreover, the rock overburden with minimum thickness of 3150 hg/cm^, 
significantly larger than for other experiments (e.g. Baksan and 1MB), re- 
duces the surface muon flux of a factor of around 10^, hence providing a 
shield to possible sources of background induced by downgoing atmospheric 
muons. These features permit a fully-automated analysis of both down- 
going and upgoing muons. The symmetric time measurement provides an 
important tool to study the detector acceptance and possible sources of 
background using the huge statistics of downward-going muons. 

The MACRO detector, described in Ref. g] is located in the HaU B of 
the Gran Sasso Laboratory. It is a large rectangular box of 76.6 x 12 x 9.3 
m^, divided longitudinally in six similar supermodules and vertically in a 
lower part (4.8 m high) and an upper part (4.5 m high), called "attico". 
The lower half of the detector is flUed with rock absorber between horizontal 
streamer tube planes, while the attico is empty and contains the electronics 
racks and work areas. There are 10 horizontal streamer planes in the 
bottom half of the detector and 4 planes on the top made of 3 cm wire 
cells and 27° stereo strip readouts. Six vertical planes of streamer tubes 
cover each lateral side of the detector. The scintillator system consists of 
around 600 tons of liquid scintillator inside 12 m long counters which make 
three horizontal layers and a vertical layer for each lateral wall. The time 
(position) resolution for muons along a scintillator box is about 500 ps 
(~ 11 cm). 

The schematic plot in Fig. illustrates the three detectable neutrino 
event topologies: Up Through, Internal Up (lU), Upgoing Stopping fis 
(UGS) and Internal Down (ID) events. The Up Through and Internal Up 
events are measured using the time-of-flight tecnique as the events cross at 
least two scintillator layers. The Upgoing Stopping /is and Internal Down 
events cross only one scintillator. The lack of time information prevents to 
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distinguish the two sub-samples. Fig. |l| shows also the parent neutrino en- 
ergy distributions for the detectable event topologies. The average neutrino 
energies are 100 GeV for Up Through niuons, 4 GeV for Upgoing Stopping 
jjLS together with Internal Down events and the same for Internal Up events. 
The Up Through jis and the Upgoing Stopping fis are externally produced 
events coming from v^j^ interactions in the rock surrounding MACRO. In 
the case of Up Through, the muons cross the detector. The absorber in 
the lower half of the detector determines an energy threshold of 1 GeV for 
vertical muons. The Internal Up and Internal Down events come from v 
interactions with vertex inside the bottom half of the apparatus. 

If the atmospheric neutrino anomaly is the result of — oscillations 
with Am^ ~ lO^'^ — 10^^ cV^ it is expected a reduction in the flux of the 
upward throughgoing muons and a distortion in the shape of their angular 
distribution. In fact, the longer pathlength of neutrinos coming from the 
nadir with respect to that coming from the horizon should cause a stronger 
reduction of the flux near the vertical than near the horizontal. On the 
other hand, for lower energy events and with these values of Am^, it is 
expected a reduction of their flux of about a factor two and no distortion 
in their angular distribution. 

The measurements of the three topologies corresponding to different 
energy ranges permit the investigation of distinct regions of the oscillation 
parameter space. Moreover, measuring two topologies with energy spectra 
in the same energy range, will allow MACRO to evaluate a double ratio, 
with consequent cancellation of most of the theoretical errors. 

2.1. Upward Throughgoing muons 

The signature for muon neutrino induced events is the versus of flight of 
muons: neutrinos can cross amounts of matter from 10 to 10^ km (the 
diameter of the Earth) and can be detected through their charged current 
interactions in the surrounding rock as upward-going muons. On the other 
hand, the neutrino induced downward-going muons cannot be discrimi- 
nated among the atmospheric muons of many orders of magnitude more 
numerous. 

The data have been collected during three periods with different detec- 
tor configurations: in the first two periods (March 1989 - November 1991 
and December 1992 - June 1993) only lower parts of MACRO were into 
acquisition; during the last period (April 1994 - November 1997) MACRO 
was put into acquisition in its final configuration including the attico. The 
analysis of the first two samples of data is described in Ref. 0. Around 3 
years of data have been analyzed with the final configuration |q] . 

In Fig. ||(a) the 1/(3 distribution is shown (1//3 is proportional to the 
measured time-of-flight) . Muons downward-going through the detector are 
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expected to have ~ +1, while muons moving upward are expected to 
have 1//3 ~ —1. This distribution is obtained by imposing several cuts to 
remove the background caused by radioactivity in narrow coincidence with 
muons, showering events and multiple muons which may result in bad time 
reconstruction. The most important requirement is that the position along 
the scintillator counter measured using the streamer tubes has to differ from 
the one from the times measured at the two ends of the counter by less 
than ±70 cm (±140 cm for nearly horizontal tracks). When a muon crosses 
three scintillator layers (almost 50% of the tracks), there is redundancy in 
the time measurement and is calculated by a linear fit of the times as a 
function of the pathlength. Tracks with poor fits are rejected. Other minor 
cuts are applied for tracks which intercept only two scintillator layers. The 
achieved rejection factor of upward throughgoing muons with respect to 
atmospheric muons is ~ 10~^. 

The requirement that throughgoing muons should cross at least 200 
g/cm^ has been applied to reduce at the percent level the background due 
to an undetected downward-going muon which interacts in the rock sur- 
rounding the detector and produces a detected upgoing particle (in most 
of the cases a pion) 0. Finally, we have observed that a large number of 
nearly horizontal upgoing muons (cosO > —0.1) come from the azimuth 
angle region between — 30°tol20°, which corresponds to a cliff in the di- 
rection of Tcramo. We exclude this angular region due to the unsufficient 
rock coverage. 

The total number of upward throughgoing muons with 1//3 in the in- 
terval [—1.25,-0.75] is 479. Based on events outside the upward-going 
muon peak, we estimate a background of 9 ± 5 events due to incorrect /3 
measurement (such as in the case of more tracks crossing the same scintil- 
lator counter) and that there are 8 ± 3 events resulting from upward-going 
charged particles produced by downward-going muons in the rock. Fi- 
nally, it is estimated that 1 1 ± 4 events are the result of f interactions in 
the bottom scintillator layer which satisfy the upward-going muon analysis 
requirements. Hence, removing the backgrounds and the internally pro- 
duced estimated events, the total number of measured upward throughgo- 
ing muons is 451. 

The expected upward-going muon flux is obtained using the Bartol u 
flux j|] , the Morfin and Tung parton distributions set Si for the calcu- 
lation of the v CC cross-sections (chosen because of the good agreement 
of the total cross section with the world average at Ei, = 100 GeV). The 
propagation of muons to the detector has been done using the energy loss 
calculation by Lohmann et al ||l^ . The total theoretical uncertainty in the 
calculation of the upward-going muon flux, obtained summing in quadra- 
ture the errors is ±17%. This error concerns the normalization of the flux, 
not the shape which does not vary changing the inputs of the calculation 
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Figure 1. On the left: diflFerent topologies induced by v interactions in or 
around MACRO: V'p Through, Internal Up (lU), Upgoing Stopping (UGS) 
jjLS, and Internal Down (ID). On the right: distributions of the parent v 
energy giving rise to the 3 different v event topologies computed by Monte 
Carlo using the same cuts applied to real data and normalizing to 1 yr of 
data taking. The average neutrino energies of the samples are: ~ 100 GeV 
for Up Through, - 4 GeV for lU and - 4 GeV for UGS and ID. 



and it is known at the level of few percents |ll) . The total number of ex- 
pected events is 612, giving a ratio of the observed to the expected events 
of 0.74 ± Om^stat ± OMQsys ± 0-Utheor- Fig. |(b) shows the zenith distri- 
bution of the measured flux of upward throughgoing muons with > 1 
GeV compared to the expected one. The dashed line represents the flux ob- 
tained in the hypothesis of t'^ oscillation with maximum mixing and 
Am^ — 0.0025 eV^. The test applied to the angular distribution gives 
= 26.1 for 8 d.o.f. (probability of 0.1% for a shape at least this different 
from the expectation) for the no-oscillation hypothesis excluding the last 
bin due to a possible contamination of atmospheric horizontal muons and 
normalizing the prediction to the data (the main theoretical error concerns 
the normalization). The best when testing the i/^ —^ Vr hypothesis 
is 15.8 for Ato^ around 0.0025 eV^ and maximum mixing (if the bound 
on the meaningful interval of sin^29 is released the best is found for 
mixing greater than 1). The independent probabilities for obtaining the 
observed number of events and the angular shape of the distribution calcu- 
lated for maximum mixing and varying Am^ are shown in Fig. |^ (plot on 
the left). The Am^ region in which the two probabilities have their max- 
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Figure 2. (a) distribution for the full detector data. The downgoing 
muon peak contains ~ 18.4 millions of events and the well separated upward 
throughgoing muon peak contains 398 events. The shaded part of the 
distribution is for the events which cross 3 scintillator layers, (b) Zenith 
distribution of the flux of upward throughgoing fj,s with i?^ > 1 GeV for 
the data and the Monte Carlo. The solid line shows the expectation for 
no-oscillations with the 17% theoretical uncertainty (shaded region). The 
dashed line is the prediction for an oscillated flux with sin^ 29 = 1 and 
Am^ = 0.0025 eV^. 



imum is similar. In Fig. ^ (plot on the right) the solid lines represent the 
iso-probability contours in the oscillation parameter space obtained com- 
bining the two independent probabilities. The curves correspond to 10% 
and 1% of the maximum probability of 17% obtained for maximum mixing 
and Am^ ~ 0.0025 eV^. The dashed lines are the contours of the confidence 
regions at the 90% and 99% c.l. based on the Monte Carlo prescription in 
Ref. jlj]. This prescription assumes that the hypothesis is correct. The 
sensitivity curve results from the preceding prescription if the data and 
Monte Carlo happened to be in perfect agreement at the maximum prob- 
ability point. The allowed regions are smaller than the sensitivity curve 
due to the fact that the absolute minimum is outside the physical re- 
gion. The same procedure has been applied considering the hypothesis of 
muon neutrino oscillations into a sterile neutrino and the maximum 
probability obtained is 2%. 
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Figure 3. On the left: probabilities for obtaining the observed MACRO 
zenith distribution shape and total number of events in the Vr oscil- 

lation hypothesis as a function of Am^ and for sin^ 29 = 1. On the right: 
probability contours from the combined probabilities in the plot on the left. 
The best probability is 17% and the iso-probability contours (solid lines) 
are shown for 10% and 1% of this value. The confidence regions (90% and 
99% c.l.) calculated according to Ref. ||l^ (dashed lines) and the sensitivity 
curve of the experiment (dotted line) are shown. 
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Figure 4. Comparison between measured and expected low energy event 
distributions vs cos9. The dashed line is obtained assuming i^^ —>■ i>r 
oscillations with parameters coming from the Up Through sample. 
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Table 1. Low energy event summary. The predictions with oscillations 



are for sin^ 26 = 1 and Ato^ = 0.0025 eV^ 



Topology 


Data 


MC (no-oscillation) 


MC(l/^ Vr) 


Up Through 


451 


612 ± lUtheor ± 37,ys 


431 ± 73ifceor ± 26,y^ 


Internal Up 


85 


144 ± 36theor ± 14s„, 


83 ± 21j(,eoi- i 8sj;s 


Int. Down + UGS 


120 


159 ± 40theor ± 16a„s 


123 it 'i^theor i 12sys 



2.2. Low energy events 

The data sample for the analyses of low energy events has been collected in 
~ 3 yr of data with the full detector. The analysis of the Internal Up events 
is similar to the Up Through analysis, with the additional requirement that 
the interaction vertex lies inside the apparatus. From a full Monte Carlo 
using the Bartol neutrino flux jsj, the v cross-sections in Rcf. and the 
Lohmann et al muon energy loss it is found that about 87% of the 
events are due to v^j^ CC interactions. The simulation has been performed 
in a large volume of rock (170 Kton) around the apparatus (5.3 Kton). 
The theoretical uncertainty estimated for the calculation is about 25%. 
The preliminary estimated systematic error on the acceptance and analysis 
requirements is ^ 10%. After background subtraction (3 events) 85 events 
are selected as Internal Up events. 

The Internal Down events and the Upgoing Stopping /is are selected 
through topological cuts: the main requirement is the presence of a track 
crossing the bottom scintillator layer at least 1 m from the detector walls. 
To reject ambiguous and/or wrongly reconstructed tracks which survived 
the software cuts, real and simulated events were randomly merged and 
visually scanned. Three different samples have been selected according to 
the choice of the minimum number of streamer plane hits. In Fig ^ we show 
the angular distribution of the 120 events (after the background subtraction 
of 5 events) which cross at least 3 streamer planes (corresponding to ~ 100 
g/cm^) and the angular distribution of the Internal Up events. 

In Tab. [| the total number of measured and expected events (no- 
oscillation and — > Vt hypothesis) are summarized for the two topologies. 
The low energy events show a uniform (as a function of cos 9) deficit of 
roughly a factor of 2 with respect to the expected ones in the no-oscillation 
hypothesis. This result is in agreement with the prediction based on the 

— > Vt oscillation hypothesis coming from the Up Through sample. 
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3. The search for dark matter 

The existence of non-baryonic dark matter has been suggested due to many 
evidences indicating that the density of matter in the universe, flM (in 
units of the critical density), is greater than 0.3. Theoretical prejudices 
for a flat universe require the density of the universe to be 11 = 1. These 
indications and hypotheses together with the upper limits on luminous 
matter {fliu7n ^ 0.01) and on baryonic mass from primordial baryonic 
nucleosynthesis (i.e. 0.005 < fife < 0.10 at 95% c.l. for 0.4 < /i < 1 |l|]) 
explain the need for non-baryonic dark matter. It is a striking coincidence 
that if a weakly interacting massive dark matter particle (WIMP) exists 
its abundance would correspond today nearly to the critical density due 
to the fact that the thermally averaged annihilation cross-section would be 
comparable to the typical weak interaction cross-section []l7t . 

If R-parity is conserved, supersymmetric (SUSY) theories would provide 
one of the most plausible candidate, the Lightest Supersymmetric Particle 
(LSP), which in most theories is the neutralino the lightest linear super- 
position of gaugino and higgsino eigenstates. In Minimal Supersymmetric 
extensions of the Standard Model (MSSM) including the GUT relation 
Ml = M2ta.T?9\Y, with 9w the Weinberg angle, the neutralino mass de- 
pends on the gaugino mass parameters Mi , on the higgsino mass parameter 
/i and on the ratio of the Higgs doublet vacuum expectation values tan [3. 
Some other parameters must be determined in order to define the processes 
induced by neutralinos, such as m^, the mass of the pseudoscalar Higgs bo- 
son, the squark and slcpton common mass mo and the trilinear parameters 
of the S'"'* family with common value A. 

The supersymmetric parameter space is constrained by collider searches 
[ p^ and lower limits on the neutralino mass are around 20-30 GeV. These 
limits are model dependent and correlated to chargino limits. In this frame- 
work, "direct" and "indirect" methods for detecting galactic halo WIMPs 
performed with underground detectors can probe complementary regions 
of the parameter space with respect to collider searches. Direct methods 
employ low-background detectors (e.g. superconductors or scintillators) to 
measure the energy deposited when a WIMP elastically scatters from a nu- 
cleus; indirect methods measure neutrinos resulting from the annihilation 
of WIMP pairs. 

WIMPs in the halo can lose energy through clastic scattering in the 
celestial bodies which they intercept and become gravitationally trapped 
inside their core when their velocity falls below the escape velocity. As 
their density increases inside the core of the body, their annihilation rate 
increases until equilibrium is achieved between capture and annihilation 
rate. High energy us are eventually produced via the hadronization and/or 
decay of the annihilation products (mostly fermion-antifermion pairs, weak 
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and Higgs bosons). Possible "traps" for WIMPs could be the Sun or the 
Earth. Many calculations concerning the flux of neutrino induced upward- 
going muons from neutralino annihilation in the Sun and the core of the 
Earth have been made |9[ |20|. Data on upward muons have been 
presented by several experiments p2| , p3[ . 



3.1. WIMP search with MACRO 

The aim of the search is to set flux limits (in absence of an evidence for 
a signal) on the upward-going muon flux from the core of the Earth and 
of the Sun. Muon flux limits are evaluated (90% c.l.) as the ratio of the 
poissonian upper limit for the number of measured events and of expected 
background events due to atmospheric i^s inside the search cone around the 
direction of the celestial body over the exposure (live time times area seen 
by the source times detection efHciency). 

For the search of WIMPs in the core of the Earth a larger sample of 517 
(collected up to Marc h '9 8) upward throughgoing muons with respect to 
the analysis in section has been used. After background subtraction the 
total resulting number of events is 487 ± 22stat while the expected events 
are 652 ± llltfi^or- As the background for this measurement overcomes the 
signal in the region around the vertical of the apparatus, we set conservative 
flux limits for WIMPs in the Earth assuming that the number of measured 
events equals the expected one [Q. The expected distribution has been 
normalized by a factor 0.85 (the ratio of measured to expected events with 
9 > 30°, which is outside the region interesting for the signal). The reason 
for this normalization is that calculations on atmospheric vs are mainly 
affected by a normalization theoretical uncertainty. The zenith angular 
distribution of measured and expected events is shown in Fig. ||(a) and the 
90% c.l. muon flux limits for an average exposure of 2620 m^ yr (the area 
decreases of about 37 m^ from 3° to 30°) are plotted in Fig. |^(c) in search 
cones increasing from 3° to 30°. Considering the results in Ref. (see 



section 2.1), the hypothesis of i/^ — > Vr with sin^ 29 = 1 and Am^ — 0.0025 
eV^ has been considered (dotted lines in Fig. ^(a) and (c)). 

As for moving sources the background is less relevant than for steady 
sources, in the case of the Sun we have used 762 events including Internal 
Up events and events which cross less than 200 g/crn^ through the ab- 
sorber. The expectation has been obtained using the local coordinates of 
the measured upward-going events and randomly chosen detection times of 
downward-going muons measured during the data taking (times are needed 
in order to evaluate equatorial coordinates). This procedure takes into ac- 
count the drifts of detection efficiency in time. Fig. ^(b) shows the dis- 
tribution of the 762 events as a function of their angular separation with 
respect to the Sun direction. In this distribution events detected during 
night fall towards -1. In Fig. ^(d) the muon flux limits are shown as a 
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Figure 5. (a) Zenith distribution of measured (circles) and expected 
(solid line) Up Through fj,. The expectation is normalized to the ratio 
data/MC for 9 > 30° equal to 0.85. (b) Distribution of measured (circles) 
and expected (solid line) upward-going /is vs the cosine of the angular 
separation from the Sun direction, (c) Muon flux limits (90% c.l.) vs the 
zenith angle. In (a) and (c) the dotted line is for in the z^^ oscillation 
hypothesis for the background of atmospheric vs with sin^ 29 — 1 and 
Am^ = 0.0025 eV^. In this case the normalization factor is 1.19. (d) Muon 
flux limits (90% c.l.) vs search angle around the direction of the Sun. 



function of 10 arbitrary search cones around the direction of the Sun for 
an exposure of ^ 890 m^ yr. Some of the limits are summarized in Tab. ^. 
The minimum muon energies have been chosen where the acceptance of 
the detector becomes independent on energy. It has an higher value for the 
Sun because tracks from the Sun direction are more slanted than the ones 
at the vertical of the detector. 
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Table 2. Selected and expected events and 90% c.l. muon flux limits for 
some of the 10 half-cones chosen around the core of the Earth (exposure 
~ 2620 m^ yr, i?^, min > 1-5 GeV) and the Sun (exposure ^ 890 m^ yr, 
E^^ min > 2 GeV) direction. 



EARTH SUN 



Cone 


Data 


Backgr. 


Flux Limit 


Data 


Backgr. 


Flux Limit 












(cm"'^ s~^) 


30° 


76 


123.6 


2.28 xlO"-^" 


56 


51.1 


6.02 xlQ-^* 


24° 


52 


81.2 


1.81 xlO-14 


33 


33.0 


3.85 xlQ-i* 


15° 


24 


32.5 


1.25 xlO"" 


11 


13.0 


2.12 xlQ-" 


9° 


10 


11.9 


7.54 X 10-15 


3 


4.6 


1.35 xlQ-i* 


3° 





1.3 


3.72 xlO-15 


2 


0.5 


1.73 xlO~i* 



3.2. Upper limits for neutralinos 

In the previous section flux limits have been calculated in arbitrary search 
cones, so that any WIMP model can be constrained by MACRO data. 
Here we assume that WIMPs are neutralinos and hence we calculate flux 
limits in the search cones which collect 90% of the signal from neutralino 
annihilation. The muons from x annihilation are distributed around the 
neutrino direction because of the ly — fj, angle due to the v CC interaction in 
the Earth and to the multiple scattering of muons along the path from the 
1/ interaction point to the detector. The v — ^ angle distribution depends 
on the V energy spectrum and hence in principle on the details of the final 
states of the annihilation process. However, the relevant parameter is Wj^. 
The branching ratios into the final states depend on the supersymmetric 
model. We estimate that the higher variation of the width of the v — ^ 
angle is of the order of 17% when an extreme model (in which only one 
channel is allowed) is used instead of a more plausible model (as the one 
used here). 

The shapes of the upward-going muon signal has been calculated using 
v fluxes from neutralino annihilation in the Sun and in the Earth given by 
Bottino et al jl^, the v CC cross sections in Ref. jlj], the muon propagation 
in Ref. ||l^ and taking into account the angular resolution of the detector. 
We have considered the x distribution in the core of the Earth |l9). The 
90% c.l. flux limits have been calculated in neutralino mass dependent 
cones which collect 90% of the signal shown in Fig. | for the Earth and 
the Sun for = 60, 100, 200, 500, 1000 GeV. In Fig. the experimental 
limits are superimposed to the fluxes of upward-going muons from Bottino 
et al |l9t calculation as a function of for a muon minimum energy of 1 
GeV (the dependence on the energy of MACRO acceptance has been taken 
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Figure 6. (a) Nadir angle distribution of fis from neutralino annihilation 
inside the Earth for different masses, (b) v — ^ angle distribution for 
neutralino annihilation in the Sun and different masses. In (a) and (b) the 
angular ranges including 90% of the signal are indicated. 

into account in the low energy region). The fluxes are calculated varying 
the model parameters (each model is represented by a dot in Fig. |^) in 
experimentally allowed ranges: 10 < |m| < 500 GcV, 10 < A/2 < 500 
GeV, 1.01 < tan/3 < 50, 65 < < 500 GeV, 150 < mo < 500 GeV, 
— 3 < A < 3, where is the mass of the pseudoscalar Higgs, mp is the 
common soft mass of all the sfermions and squarks, A is the common value 
of the trilinear coupling in the superpotential for the third family, with 
central values of allowed interval for cosmological parameters (rms velocity 
of X in the halo — 270 km s~^, x escape velocity in the halo = 650 km s~^, 
velocity of the Sun around the galactic center — 232 km s~^, local dark 
matter density pioc = 0.5 GeV cm"'^ and minimal value for rescaling the 
neutralino relic abundance flhf^^^ — 0.03). In fact, Bottino et al. assume: 

if ilft-l > {flh'^)min Px = Ploc 

if nh\ < {nh'^)min Px = Ploc X f7ft.|/(f]/l2)„„ . 

where (Q,h?')min — 0.03. The fluxes lying above the experimental limit are 
ruled out as possible SUSY models by this measurement. Fig. |^ indicates 
that indirect search can have good prospects particularly at high x masses. 

4. Conclusions 

The upward throughgoing muon MACRO data set is in favor of a Vr 
hypothesis with Am^ ~ 0.0025 eV^ and maximum mixing with maxi- 

13 



mum probability of 17% against the 0.1% probability obtained for the 
no-oscillation hypothesis. These oscillation parameters are very close to 
the ones preferred by the SuperKamiokande experiment [ p5| . The hypoth- 
esis found is in agreement with the first results concerning the MACRO 
low energy neutrino events. However, the study of the shape of the upward 
throughgoing muon angular distribution gives a maximum probability of 
only 4.6%. This could be due to a statistical fluctuation or to some hidden 
physics. Having performed numerous systematic checks, we exclude any 
effect on the anomalous zenith distribution due to the detector. 

A search for a WIMP signal has been performed using throughgoing 
upward /iS coming from the Earth core and from the Sun. At this time we 
have observed no signal for this signature and we have calculated upper 
limits for the flux of upward /is from WIMPs. These flux limits have been 
compared with predictions of SUSY models. Our data have begun to rule 
out significant portions of the parameter space for WIMPs annihilating in 
the Earth and in the Sun and currently limits for the Earth are the most 
stringent limits of all those from "indirect" experiments j2^, |2|] . 
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